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have  also  been  obtained. 
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Included  the  study  of  linear  empirical  Bayes  estimators,  prediction  Intervals 
for  the  Inverse  Gaussian  distribution,  nonparametrlc  hazard  rate  estimation 
under  censoring,  nonparametrlc  Inference  for  step-stress  accelerated  life  tests 
under  censoring,  discrete  failure  models,  simultaneous  confidence  Intervals 
for  pairwise  differences  of  normal  means,  and  optimal  designs  for  comparing 
treatments  with  a  control.  In  addition,  reliability  estimation  when  the  cause 
of  failure  may  be  only  partially  known  and  problems  giving  rise  to  Gompertzlan 
failure  models  were  studied. 
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1.  SUMMARY 


This  document  reports  the  work  performed  eind  other  research  activities  of 
the  investigators  during  the  funding  period  from  June  1,  1984,  to  May  31,  1987, 
xinder  grant  number  NIPR  ARO  139-85. 

^  The  main  objectives  of  this  research  have  been  the  development  of  smooth 
nonparametric  estimators  of  quzmtile  functions  frcsn  right-censored  data  and  the 
further  study  of  smooth  density  estimators  from  censored  observations.  In 
particular,  kernel-type  and  “generalized”  quantile  estimators  have  been 
obtained  under  censoring  v^ich  give  better  estimates  of  percentiles  of  the 
lifetime  distribution  than  the  usual  product-limit  quantile  estimator. 
Asymptotic  properties  of  these  quantile  estimators  have  been  developed, 
including  asymptotic  normality  and  consistency.  Also,  the  me£ui-square 
convergence  and  mean  squared  error  of  the  kernel  quantile  estimators  have  been 
investigated.  A  data-based  procedure  for  selecting  the  bandwidth  has  been 
studied  using  the  bootstrap,  and  approximate  confidence  intervals  for  the  true 
quantile  have  been  proposed  using  bootstrap  estimates  of  the  sampling 
distribution.  Theoretical  results  on  the  optimal  bandwidth  selection  for 
kernel  density  estimators  under  random  right  censorship  have  also  been 
obtained.  — 

New  results  in  several  other  problem  areas  were  also  developed.  These 
inclixled  the  study  of  linear  empirical  Bayes  estimators,  prediction  intervals 
for  the  inverse  Gaussian  distribution,  nonparametric  hazard  rate  estimation 
under  censoring,  nonparametric  inference  for  step-stress  accelerated  life  tests 
under  censoring,  discrete  failure  models,  reliability  estimation  when  cause  of 
failure  is  partially  known,  Gompertzieui  failure  models,  simultaneous  confidence 
intervals  for  pairwise  differences  of  normal  means,  and  optimal  designs  for 
comparing  treatments  with  a  control.  All  of  these  results  should  prove 
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£ruitful  in  the  assessment  of  reliability  and  maintenance  policies  of  various 
types  of  equipment. 


2.  RESEARCH  CBJECTIVES 

Ihe  research  objectives  of  this  project  can  be  divided  into  the  following 
categories: 

A.  Nonparametric  estimation  of  queuitiles  under  censoring; 

B.  Kernel  density  estimation  frcmt  randomly  right-censored  data; 

C.  Nonparametric  estimation  from  step-stress  accelerated  life  tests; 

D.  Linear  enpirical  Bayes  estimation  of  mean  time  to  failure;  and 

E.  Miscellaneous  problems  in  reliability. 

Ihe  specific  research  problems  that  were  considered  during  this  project  period 
in  these  areas  will  be  briefly  outlined  below.  The  results  obtained  will  be 
described  in  Section  3  of  this  report. 

A.  Nonparametric  estimation  of  quantiles  under  censoring.  One  of  the 
main  objectives  of  this  research  was  to  develop  smooth  noi^rametric  estimates 
of  the  qu^ultile  fvinction  (or  percentiles)  from  randomly  right-censored  data. 

In  particular,  the  kernel-type  quantile  estimators  were  to  be  studied  with 
respect  to  small-sanple,  as  well  as  asymptotic,  properties.  Also,  a 
modification  of  the  kernel  estimator  which  allowed  the  bandwidth  to  be  a 
function  of  the  censored  data  was  to  be  studied.  A  different  approach  to 
smooth  quantile  estimation  was  taken  in  developing  a  "generalized"  product- 
limit  quantile  estimator.  This  estimator  is  essentially  an  average  of  product- 
limit  qu£uitiles  frcmt  all  possible  svbsamples  of  a  given  size.  These  types  of 
estimators  smooth  the  product-limit  queuntile  function  Q^(p)  * 
inf{t:  denotes  the  product-limit  estimator  of  the 

lifetime  distribution  F^,  and  are  better  than  Q^(p)  in  the  sense  of  smaller 


-3- 


meaun  squared  errors  in  many  situations.  Thus,  such  estimators  can  provide  more 
accurate  estimates  of  percentiles,  aiding  in  the  assessment  of  equipment 
reliability. 

B.  Kernel  density  estimation  from  randomly  right-censored  data.  Another 
objective  of  this  project  during  the  reporting  period  was  to  determine  a  method 
of  "optimal"  bandwidth  selection  for  kernel-type  density  estimators  under 
random  right-censorship,  with  respect  to  scxne  optimality  criterion.  Itie  kernel 
density  estimators  of  Blum  and  Susarla  (Multivariate  Analysis-V,  1980)  and 
Foldes,  Rejto  and  Winter  (Periodica  Mathematica  Hungarica,  1980)  were  studied. 

C.  Nonparametric  estimation  from  step-stress  accelerated  life  tests. 
Usually  in  accelerated  life  tests,  items  are  put  on  test  at  several  accelerated 
stresses  and  observed  until  failure  or  until  the  test  terminates.  A  more 
realistic  method  is  to  place  items  on  test  and  allow  the  stress  levels  to 
change  at  preassigned  times.  This  will  result  in  more  failures  in  general. 

The  objective  of  the  research  in  this  area  was  to  obtain  a  strongly  consistent 
nonparametric  estimator  of  the  unknown  life  distrilsution  at  the  normal  use 
stress.  A  rescaling  technique  was  utilized  to  study  this  problem. 

D.  Linear  empirical  Bayes  estimation  of  mean  time  to  failure.  In 
estimating  the  mean  lifetimes  of  n  independent  problems,  it  is  hoped  that  the 
linear  empirical  Bayes  estimators  will  converge  to  the  optimal  linear  Bayes 
estimators  with  reasonable  rapidity.  This  convergence  was  to  be  studied  here. 
Also,  another  objective  in  this  area  was  to  identify  a  natural  class  of 
estimators  between  the  linear  Bayes  estimators  and  the  Bayes  estimators. 

E.  Miscellaneous.  Problems  in  several  other  areas  of  investigation 
related  to  the  general  problem  of  reliability  assessment  were  to  be  studied. 
These  inclixle  hazard  rate  estimation  under  censoring,  prediction  intervals, 
discrete  failure  models,  simultaneous  inference,  euid  optimal  design  for 
treatment  comparisons  with  a  control. 
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3.  RESEARCH  ACCOMPLISHMENTS 

In  this  section  a  sxibstantive  statement  of  the  progress  and  significant 
accomplishments  towards  achieving  the  research  objectives  for  this  project 
outlined  in  section  2  will  be  given.  The  specific  research  papers  referred  to 
will  be  listed  in  section  4  of  this  report,  v^ich  shows  the  cumulative  results 
of  the  research  project. 

A.  Nonparametric  estimation  of  quantiles  under  censoring.  Based  on 
right-censored  cata  from  a  lifetime  distribution  F^,  a  kernel-type  estimator  of 
the  quemtile  function  Q°(p)  -  inf{t:  ^^(t)  >  p},  0  <  p  <  1,  was  proposed  by 
Padgett  in  paper  number  [5]  listed  in  section  4.  The  estimator  is  defined  by 
Qn(p)  ■  h~^jQQ^(t)K( (t-p)/hj^)dt,  vdiich  is  smoother  than  the  usual  product-limit 
function  Qj^(p)  *  inf{t:  >  p},  where  denotes  the  product-limit 

estimator  of  F^.  Under  the  random  censorship  model  auid  general  (but 
nonrestrictive)  conditions  on  h^,  K,  and  F^,  it  was  shown  in  [5]  that  Qj^(p)  is 
strongly  consistent.  In  addition,  it  was  shown  that  an  approximation,  Q^,  to 
is  asynptotically  equivalent  to  it  with  probability  one.  The  general 
conditions  for  consistency  included  that  h^  -»  0  as  n  -♦  •,  that  K  be  a  symmetric 
probedaility  density  with  finite  support,  and  that  F^  satisfy  certain  continuity 
conditions. 

During  this  project,  the  estimator  Qj^(p)  was  studied  further  and  in  much 
more  detail.  In  paper  [7]  it  was  shown  that,  under  conditions  similar  to  those 
required  for  strong  consistency,  n^(Q^(p)  -  Q°(p))  is  asynptotically  normally 

L 

distributed,  provided  n^h^  -»  0.  This  condition  on  h^  cem  be  replaced  by 
n^h^  -  0  by  using  a  slightly  different  proof.  By  centering  with  the  queuitity 
Q°(p,h)  -  h”^jQQ°(t)K( (t-p)/hj^)dt  instead  of  Q°(p),  an  asymptotic  normality 
result  can  be  obtained  without  the  above  condition  on  the  rate  of  convergence 
of  h^  to  zero.  Also,  using  this  type  of  centering,  the  asymptotic  normality  of 


the  approximation  was  shown,  ihese  last  two  results  are  contained  in  paper 
number  [9].  In  addition,  in  (9],  mean  convergence  of  Qj^(p)  and  Q^(p)  were 
obtained  and  the  mean  square  convergence  was  shown.  In  particular,  an  upper 
bound  on  the  rate  of  mean  square  convergence  of  Qj^(p)  to  Q°(p)  was  shown  to  be 

_u  -5/4  3/4  -i*  2 

g(n,h^)  -  o(hj^  +  hj^^n  '  +  (log  n/n)  ^  +  n  c(l-p)  ),  vrtiere  c  is  a  constant. 

An  example  of  a  bandwidth  sequence  satisfying  g(n,h^)  -»  0  as  n  ->  *  is 
h_  =  c  n~^  for  0  <  S  <  5/2  where  c_  >  0  is  a  bounded  sequence.  The  value  of  h 
minimizing  g(n,h)  was  found  to  be  of  the  form  h^  -  Cn~^^^.  This  rate  is  not 
sharp,  but  seems  to  be  close  to  optimal  for  large  n. 

New  asymptotic  expressions  for  the  meem  squared  errors  of  the  product-limit 
quantile  estimator  and  the  kernel-type  quantile  estimator  were  presented  in 
paper  [14].  From  these  expressions  a  conparison  of  these  two  quemtile 
estimators  with  respect  to  their  mean  squared  errors  was  given.  In  particular, 
the  comparison  showed  that  for  sufficiently  large  n  and  small  enough  h,  the 
kernel  estimator's  mean  squared  error  was  less. 

Since  the  exact  small-sanple  mean  squared  error  has  not  been  calculated,  an 
expression  for  "optimal"  h^  minimizing  the  roe2un  sqxaared  error  is  not  availedale. 
Hence,  a  study  of  the  effect  of  varying  the  beuidwidth  was  done  by  Monte  Carlo 
simulations  in  paper  [10].  Five  ccamnon  life  distributions,  two  different 
censoring  distributions,  and  three  different  kernel  functions  were  used  in  the 
extensive  sumulations .  The  results  indicate  that  at  fixed  p  the  kernel 
estimator  Qjj(p)  has  smaller  mean  squared  error  than  the  product-limit  quantile 

A 

estimator  Qj^(p)  for  a  range  of  values  of  the  bandwidth  h^.  In  addition,  in 
[10]  a  method  of  selecting  em  "optimal"  beuidwidth  value,  in  the  sense  of  small 
estimated  mean  squared  error,  based  on  the  bootstrap,  was  investigated.  The 
results  were  consistent  with  the  simulation  study.  The  bootstrap  was  also  used 
to  obtain  interval  estimates  of  0°(p)  after  determining  a  value  of  h^  to  use  in 
computing  Q„(p)» 


-6- 


In  paper  [19]  a  mcxiification  o£  the  kernel  qu2uitile  estimator  was  proposed 
in  vdiich  the  bandwidth  h  was  a  function  of  the  censored  data.  The  adveu'tage  of 
this  estimator  was  that  the  data  play  a  role  in  the  degree  of  smoothing  while 
retaining  the  desirable  features  of  the  kernel  estimator.  Consistency  of  the 
estimator  was  proven,  and  it  was  shown  that  a  nearest-neighbor-type  selection 
of  the  beuidwidth  gave  a  consistent  estimator. 

In  [20],  the  smooth  quantile  estimation  problem  from  censored  data  was 
approached  by  integrating  a  kernel  density  estimate  to  obtain  a  smooth 
estimate,  Fq,  of  the  lifetime  distribution.  Ihen  a  solution  Xp  from  FQ(Xp)  -  p 
was  obtained  as  an  estimate  of  the  p^  quantile.  Asymptotic  properties  of  Xp 
were  shown  and  some  computational  results  were  given. 

Finally,  in  paper  [17],  for  right-censored  data  a  smooth  alternative  to  the 
product-limit  estimator  as  a  nonparametric  queuitile  estimator  was  proposed. 

Ihis  estimator,  called  a  "generalized  product-limit  quantile  estimator,"  was 
obtained  by  averaging  appropriate  subsample  product-limit  quantiles  over  all 
subsanples  of  a  fixed  size.  Under  the  random  censorship  model  and  some 
conditions  on  the  lifetime  distribution,  it  was  shown  that  the  estimator  was 
consistent  and  had  the  same  asymptotic  normal  distribution  as  the  product-limit 
quekntile  estimator.  A  small  Monte  Carlo  simulation  study  showed  that  some 
values  of  the  subsanple  size  existed  for  which  this  "generalized"  estimator 
performed  better  them  the  usual  product-limit  quemtile  estimator  in  the  sense 
of  estimated  (from  the  simulations)  meem  squared  errors. 

B.  Kernel  density  estimation  from  randomly  right-censored  data.  The 
problem  of  asynptotically  optimal  selection  of  the  bandwidth  sequence  for 
kernel  density  estimators  under  censoring  was  studied  in  paper  number  [13]. 

—  1  CD  ^ 

The  usual  kernel  density  estimator  -  h^^  /_cbK(  (x-t)/hj^)dFj^(t)  and  the 

Blum-Susarla  estimator  given  by  f*(x)  -  [nhj^H*(x)  ]”^  ^j«l*^^^*~^j^'^n^*^[Aj-l] 


were  considered,  vdiere  H*  denotes  an  estimator  of  the  censoring  survival 
function,  (Xj,Aj),  j-l,...,n,  denote  the  right-censored  data,  and  denotes 
the  indicator  function  of  the  event  [ • ) .  Two  inportant  contributions  to  the 
theory  of  bandwidth  selection  for  kernel  density  estimators  under  right 
censorship  were  made  in  this  paper.  First,  an  asymptotic  representation  of  the 
integrated  squared  error  of  each  estimator  in  terms  of  easily  understood 
variance  emd  squared  bias  conponents  is  given.  Second,  it  is  shown  that  if  the 
bandwidth  is  chosen  by  the  data-based  method  of  least  squares  cross-validation, 
then  it  is  asynptotically  optimal  in  a  compelling  sense.  The  optimal  bandwidth 
is  shown  to  be  ^ulalogous  to  that  for  the  case  of  no  censoring.  A  by-product  of 
the  first  part  is  an  interesting  comparison  of  the  two  kernel  density 
estimators  vdiich  shows  that  the  difference  between  the  two  is  typically  not 
negligible. 

C.  Nonparametric  estimation  from  step-stress  accelerated  life  tests.  In 
accelerated  life  testing,  items  are  subjected  to  greater  stress  than  that  under 
the  .normal  use  conditions  and,  from  the  resulting  failure  data,  an  estimate  of 
the  lifetime  distribution  (or  other  inferences)  under  the  normal  use  conditions 
is  obtained.  Most  of  the  published  work  on  inference  from  accelerated  life 
tests  require  at  least  one  of  the  following  restrictions: 

(i)  that  every  item  which  is  subjected  to  an  overstress  is  observed  until  it 
fails  or  is  removed  (censored)  under  constant  application  of  the  stress,  or 

(ii)  that  the  lifetime  distribution  of  the  items  under  every  accelerated  stress 
level  is  assumed  to  be  known  except  for  the  underlying  pareuneters.  A  more 
realistic  approach  allows  the  stress  on  an  unfailed  item  to  change  at  a 
preassigned  test  time.  This  approach  is  used  in  government  defense  studies  and 
is  referred  to  as  step-stress  testing. 
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In  paper  number  [15]  step-stress  testing  is  considered  from  a  nonparametric 
•  point  of  view  vrtien  the  observations  have  been  arbitrarily  right  censored, 

>  relaxing  both  restrictions  (i)  and  (ii)  eibove.  Thus,  no  assumptions  are  made 

about  the  particular  form  of  the  underlying  lifetime  distribution.  Using  a 
rescaling  method,  an  estimator  of  the  life  distribution  under  the  normal  use 

I 

conditions  is  obtained.  The  estimator  is  shown  to  be  strongly  consistent. 

D.  Linear  empirical  Bayes  estimation  of  mean  time  to  failure.  In  paper 
[11],  it  has  been  successfully  established  that  not  only  the  linear  empirical 

I 

I  Bayes  estimator  approaches  the  optimal  linear  estimator  with  probability  one, 

but  also  in  the  quadratic  meeui  with  a  rate  of  convergence  inversely 
proportional  to  the  nxjmber  of  problems  considered.  The  same  is  true  for  the 
meeui  squared  error  of  the  estimators.  Indeed,  the  limiting  constant  has  been 
explicitly  evaluated,  v^ich  will  provide  a  guideline  for  the  consideration  of  a 
moderate  nximber  of  problems.  Effort  will  be  continued  to  study  the  moderate 
sample  properties  of  the  linear  empirical  Bayes  estimators,  and  to  identify  a 
natural  and  appropriate  intermediate  class  of  estimators. 

Progress  has  also  been  made  on  an  investigation  into  identification  of  a 
natural  and  appropriate  intermediate  class  of  estimators.  Specifically,  the 
class  of  quadratic  estimators  has  been  singled  out  and  some  preliminary 
properties  have  been  shown.  Further  work  needs  to  be  done  to  obtain 
comprehensive  and  definitive  results. 

E.  Miscellaneous.  In  other  work,  an  invited  paper  [4]  was  written  on 
nonparametric  methods  for  hazard  rate  estimation  from  censored  data.  This 
paper  first  reviewed  the  available  methods  in  the  literature  and  then  presented 
some  new  results  on  a  kernel-type  estimator  of  the  hazard  rate  under  remdom 
right  censoring.  Ihe  Koziol-Green  (or  proportional  hazards)  model  was  assumed 
for  asymptotic  results  on  the  bias  and  mean  squared  error  of  the  kernel 
estimator. 
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In  paper  number  [8],  a  conparison  was  given  of  two  prediction  intervals  for 
a  future  observation  from  the  inverse  Gaussian  distribution.  This  distribution 
arises  as  a  lifetime  model  vdiich  is  an  alternative  to  the  lognormal 
distribution.  One  prediction  interval  is  exact  amd  the  other  is  an  approximate 
one.  Surprisingly,  the  conparison  showed  that  the  approximate  interval  is 
superior  in  general,  since  it  always  yields  a  two-sided  interval  and  has  high 
coverage  probability  as  well  as  small  expected  length. 

In  paper  number  [16]  a  regenerative  sampling  plan  is  developed  for  the 
sequential  conparison  of  two  populations  having  positive  integral  response,  it 
is  designed  to  be  both  em  extension  euid  an  improvement  of  the  play-the-winner 
rules  for  binary  trials  in  the  sense  that  a  much  wider  variety  of  responses  is 
allowed,  the  fraction  of  inferior  selections  approaches  zero,  eind  the  play-the- 
winner  rule  is  contained  as  a  special  case.  Almost  sure  convergence  and  moment 
convergence  in  the  p^  order  is  studied  for  the  fraction  of  inferior  selections 
emd  for  a  meucimum  likelihood  estimator  of  the  mean  response.  A  conditional 
test  of  hypothesis  is  given  for  the  binary  case.  These  results  have 
applications  in  clinical  trials  and  in  quality  control  as  well  as  in  a  variety 
of  other  sequential  situations. 

In  paper  number  [18],  a  simple  procedure  is  proposed  to  determine  a  sample 
size  for  estimating  the  roeeui  weight  of  items  in  a  problem  of  obtaining  a  large 
number  of  items.  The  situation  is  as  follows:  It  is  desired  to  obtain  a  large 
number,  N  ,  of  items  for  which  individual  counting  is  impractical,  but  one  can 
demand  a  batch  to  weigh  at  least  w  units  emd  hope  that  the  number  of  items  in 
the  batch  is  close  to  the  desired  number  N  .  If  the  items  have  mean  weight  0, 
it  is  reasonable  to  have  w  equal  to  ew  when  0  is  known.  When  0  is  unknown, 
one  can  take  a  sample  of  size  n,  not  bigger  than  N^,  to  estimate  0  by  a  good 
estimator  0^  and  set  w  equal  to  The  proposed  procedure  determines  the 
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sanple  size  to  be  the  integer  closest  to  pCN  ,  v^ere  C  is  a  function  of  the 

cost  coefficients  if  the  coefficient  of  variation  p  is  known.  It  is  known  to 

be  optimal  in  some  sense.  If  p  is  unknown,  a  simple  sequential  procedure  is 

proposed  for  vdiich  the  average  saitple  number  is  shown  to  be  asynf)totically 

equal  to  the  optimal  fixed  sample  size.  When  the  weights  are  assumed  to  have  a 

gamma  distribution  given  0  and  0  has  a  prior  inverted  gamma  distribution,  the 

optimal  sample  size  in  some  sense  Ccin  be  found  to  be  the  nonnegative  integer 

2 

closest  to  pCN  +  p  A(pC-l),  vdiere  A  is  a  known  constant  given  in  the  prior 
distribution. 

In  paper  [21],  it  was  shown  that  the  Gompertzian  or  Weibull  (with 
increasing  failure  rate,  IFR)  distribution  cannot  arise  from  any  mixture  of 
exponential  distributions.  The  proof  of  this  used  well-known  results  in 
reliedaility  theory  (but  not  well-known  in  the  biology  of  aging).  Therefore,  if 
the  population  is  Weibull  (IFR)  or  Gonpertzian,  has  bathtub-shaped  failure 
rate,  or  has  a  failure  rate  that  increases  strictly  over  scxne  interval,  then 
there  must  be  units  that  individually  do  not  have  exponential  life 
distribution.  This  result  is  inportant  in  the  analysis  of  reliedaility  data. 


CUMULATIVE  LIST  OF  WRITTEN  PUBLICATiaaS 


In  this  section,  the  research  papers  that  have  been  written  under  this 
grant  are  listed.  They  are  divided  into  three  categories:  In  print,  accepted 
for  publication,  and  submitted.  Copies  of  the  manuscripts  and  reprints  have 
been  forwarded  to  the  Program  Manager  as  they  were  submitted  or  published. 

A.  In  Print 


[1]  Andre  M.  Lubecke  and  W.  J.  Padgett,  Nonparametric  maximum  penalized 
likelihood  estimation  of  a  density  from  arbitrarily  right-censored 
observations.  Communications  in  Statistics-Theory  &  Methods  14 
(1985),  257-271. 


;21  J.  D.  Spurrier  and  S.  P.  Isham,  Exact  simultaneous  confidence 

intervals  for  pairwise  comparisons  of  three  normal  means.  Journal  of 
the  Americeun  Statistical  Association  80  (1985),  438-442. 
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[3]  W.  J.  Padgett  and  J.  D.  Spurrier,  On  discrete  failure  mcxlels.  IEEE 
Transactions  on  Reliadbility  R-34  (1985),  253-256. 

[4]  D.  T.  McNichols  and  W.  J.  Padgett,  Nonpar ametric  methods  for  hazard 
rate  estimation  from  right-censored  sanples.  Journal  of  the  Chinese 
Statistical  Association  23  (1985),  1-15. 

[5]  W.  J.  Padgett,  A  kernel-type  estimator  of  a  quemtile  function  from 
right-censored  data.  Journal  of  the  American  Statistical  Association 
^  (1986),  215-222. 

[6]  J.  D.  Spurrier  and  D.  Edwards,  An  asymptotically  optimal  subclass  of 
balemced  treatment  incomplete  block  designs  for  conparisons  with  a 
control.  Bioroetrika  73  (1986),  191-199. 

[7]  Y.  L.  Lio,  W.  J.  Padgett  and  K.  F.  Yu.  On  the  asymptotic  properties 
of  a  kernel-type  quemtile  estimator  from  censored  samples.  Journal  of 
Statistical  Planning  and  inference  (1986),  169-177. 

[8]  W.  J.  Padgett  and  S.  H.  Tsoi,  On  prediction  intervals  for  future 
observations  from  the  inverse  Gaussian  distribution.  IEEE 
Transactions  on  Reliability  R-35  (1986),  406-408. 

[9]  Y.  L.  Lio  and  W.  J.  Padgett,  Some  convergence  results  for  kernel-type 
quantile  estimators  under  censoring.  Statistics  and  Probetbility 
Letters  5  (1987),  5-14. 

110)  W.  J.  Padgett  and  L.  A.  ihonibs.  Smooth  nonparametric  quantile 
estimation  under  censoring:  Simulations  cuid  bootstrap  methods. 
Communications  in  Statistics-Simulation  and  Computation  15  (1986), 
1003-1025. - 

(11]  K.  F.  Yu,  On  the  bounded  regret  of  empirical  Bayes  estimators. 

Communications  in  Statistics-Theory  and  Methods  ^  (1986),  2391-2403. 


B.  Accepted  for  Publication 

[12]  W.  J.  Padgett,  Nonparametric  estimation  of  density  and  hazard  rate 
functions  vdien  samples  are  censored.  Handlaook  of  Statistics,  Vol.  7, 
P.  R.  Ktishnaiah,  Ed.  (to  appear). 

[13]  J.  S.  Marron  and  W.  J.  Padgett,  Asymptotically  optimal  bandwidth 
selection  for  kernel  density  estimators  from  randomly  right-censored 
samples.  Annals  of  Statistics  (to  ajpear,  Dec.  1987). 

[14]  Y.  L.  Lio  and  W.  J.  Padgett,  On  the  wean  squared  error  of 
nonparametric  quantile  estimators  under  random  right-censorship. 
Communications  in  Statistics,  Theory  and  Methods  (to  appear). 


C.  Submitted  for  Publication 

[15]  D.  T.  McNichols  and  W.  J.  Padgett,  Inference  for  step-stress 

accelerated  life  tests  under  arbitrary  right-censorship.  Being 
revised  for  Journal  of  Statistical  Planning  and  Inference. 


[16]  S.  D.  Durh2un  emd  K.  F.  Yu,  Regenerative  sanpling  and  monotonic 
branching  processes.  Submitted  to  Annals  of  Statistics. 

[17]  Y.  L.  Lio  and  W.  J.  Padgett,  A  generalized  quantile  estimator  under 
censoring.  Submitted  to  Journal  of  the  American  Statistical 
Association. 

[18]  K.  F.  Yu,  On  determining  the  weight  for  obtaining  a  large  number  of 
items.  Submitted  to  Journal  of  the  American  Statistical  Association. 

[19]  Y.  L.  Lio  auid  W.  J.  Padgett,  A  modified  kernel  quzuitile  estimator 
under  censoring.  Submitted  to  Canadian  Journal  of  Statistics. 

[20]  W.  J.  Padgett  and  L.  A.  Thombs,  A  smooth  nonparametric  queuitile 
estimator  from  right-censored  data.  Submitted  to  Techncanetrics. 

[21]  F.  M.  Guess  and  M.  Witten,  A  population  of  exponentially  distributed 
individual  lifespeins  ceinnot  lead  to  Gonpertzieui  or  to  Weibull  (with 
increasing  mortality  rate)  dynamics.  Submitted  to  Science . 

5.  PROFESSIONAL  PERSONNEL  ASSOCIATED  WITH  THE  RESEARCH  EFFORT 

In  addition  to  the  principal  investigator,  W.  J.  Padgett,  several  co¬ 
investigators  have  been  partially  supported  by  this  grant  since  June  1,  1984. 
John  D.  Spurrier  was  partially  supported  during  the  period  from  June  1,  1984, 
to  May  31,  1985,  and  Kai  F.  Yu  was  partially  supported  from  June  1,  1985,  to 
May  31,  1987.  During  the  summer  of  1985,  Diane  T.  McNichols,  Department  of 
Statistics,  Virginia  Polytechnic  Institute  and  State  U.,  was  partially 
supported.  Frank  M.  Guess  was  partially  supported  as  a  co-investigator  from 
July  1,  1986,  to  May  31,  1987.  Also,  two  graduate  students  have  been 
supported,  Andre  M.  Lubecke  during  the  summer  of  1984,  eund  Y.  L.  Lio  from 
August  16,  1984,  to  May  31,  1987.  Ms.  Lubecke  completed  the  Ph.D.  degree  in 
August,  1985,  and  Mr.  Lio  will  obtain  the  Ph.D.  degree  in  August,  1987. 


6.  INTERACTIONS 

The  investigators  attended  several  meetings  and  conferences  and  gave 

(invited  and  contributed)  talks  as  follows; 

i)  W.  J.  Padgett,  J.  D.  Spurrier  attended  the  Conference  on  Reliability 
and  Quality  Control,  University  of  Missouri,  Columbia,  MO,  June  4-8, 
1984. 
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ii)  W.  J.  Padgett  (joint  with  A.  M.  Lubecke),  "Nonparametric  maximum 

penalized  likelihood  estimation  of  a  density  from  arbitrarily  right- 
censored  data."  Annual  Meeting  of  IMS  in  Tahoe  City,  CA,  August  21- 
24,  1984. 

iii)  J.  D.  Spurrier,  "The  training  of  statistical  consultants  at  the 

University  of  South  Carolina"  (invited).  Joint  Statistical  Meetings, 
Philadelphia,  PA,  August  1984. 

iv)  D.  Edwards,  J.  D.  Spurrier,  "An  asymptotically  optimal  subclass  of 
BTIB  designs  for  conparisons  with  a  control."  Joint  Statistical 
Meetings,  Philadelphia,  PA,  August,  1984. 

v)  J.  D.  Spurrier  and  W.  J.  Padgett,  "On  discrete  failure  models." 

Spring  Regional  Meeting  of  the  Biometric  Society  (ENAR)  and  ASA  in 
Raleigh,  NC,  March  25-27,  1985. 

vi)  W.  J.  Padgett  was  the  invited  discusseuit  for  the  Session  on 

ReliaJsility  and  Life  Testing,  Spring  Regional  Meeting  of  the  Biometric 
Society  (ENAR)  and  ASA  in  Raleigh,  NC,  March  25-27,  1985. 

vii)  W.  J.  Padgett,  "Some  estimation  procedures  for  queuitile  emd  density 
functions  from  censored  data."  AFOSR  Workshop  on  Reliability,  Luray, 
VA,  May  29-31,  1985. 

viii)  W.  J.  Padgett  gave  an  invited  talk  on  "Nonparametric  density  and 

failure  rate  estimators  vhen  samples  are  censored"  at  the  SREB-ASA 
Summer  Research  Conference  on  Statistics,  Boone,  NC,  June  16-21,  1985. 

ix)  W.  J.  Padgett,  "A  kernel  type  estimator  of  a  quantile  fvinction  from 
right-censored  data."  Joint  Statistical  Meetings,  Las  Vegas,  NV, 
August  5-8,  1985. 

x)  K.  F.  Yu,  "A  note  on  the  bounded  regret  of  empirical  Bayes 

estimators,"  191st  Meeting  of  the  IMS,  Stoney  Brook,  NY,  June  12-13, 
1985. 

xi)  K.  F.  Yu,  "On  the  bounded  regret  of  empirical  Bayes  estimators"  Joint 
Statistical  Meetings,  Las  Vegas,  NV,  August  5-8,  1985. 

xii)  K.  F.  Yu,  "On  a  renewal  theorem."  Centenary  Session  of  the  ISI, 
Amsterdam,  The  Netherlands,  August,  1985. 

xiii)  K.  F.  Yu  organized  and  chaired  a  session  on  Adaptive  Statistics  at  the 
Synposium  on  Adaptive  Statistical  Inference  emd  Related  Topics, 
Brookhaven  National  Laboratory,  NY,  June  8-11,  1985. 

xiv)  K.  F.  Yu  attended  the  NSF-CBMS  Conference  on  Nonparametric  Priors, 
university  Park,  PA,  June  3-7,  1985. 

xv)  W.  J.  Padgett  (with  Y.  L.  Lio),  "Some  convergence  results  for  kernel- 
type  quantile  estimators  from  censored  data,"  Spring  Statistics 
Meetings  (ASA/ENAR),  Atlanta,  GA,  March  16-19,  1986. 

xvi)  W.  J.  Padgett  gave  invited  talk  on  "Some  results  for  kernel-type 
quantile  estimators  from  right-censored  data"  at  the  U.S.  Army 


Research  Office  Workshop  on  Life  Testing,  Washington,  DC,  March  9-12, 
1986. 


i)  W.  J.  Padgett  and  K.  F.  Yu  attended  the  8th  Anniaal  Meeting  of  the  S.C. 
Chapter  of  ASA,  Columbia,  SC,  April  25,  1556. 

i)  K.  F.  Yu  attended  the  Spring  Statistics  Meetings,  Atlanta,  GA,  March 
16-19,  1986. 

v)  W.  J.  Padgett  gave  an  invited  talk  on  "Asymptotic  results  for  kernel 
quantile  estimators  under  censoring"  at  the  Charlotte  Mathematics 
Conference  (AMS/MAA),  Charlotte,  NC,  on  October  17-18,  1986. 

x)  W.  J.  Padgett  and  L.  A.  Thombs  presented  a  talk  at  the  ASA  Winter 
Conference  on  "Smooth  nonparametric  queuitile  estimation  under 
censoring:  Simulations  and  bootstrap  methods"  in  Orlando,  FL,  Jamuary 
7-9,  1987. 

;i)  Y.  L.  Lio  amd  W.  J.  Padgett  gave  a  contributed  paper  on  "A  generalized 
quaintile  estimator  under  censoring"  at  the  IMS  Eastern  Regional 
Meeting,  Blacksburg,  VA,  May  27-29,  1987. 

i)  W.  J.  Padgett  and  K.  Yu  attended  the  Conference  in  Honor  of  I.  J.  Good 
on  May  25-26,  1987,  and  the  IMS  Eastern  Regional  Meeting,  May  27-29, 
1987,  in  Blacksburg,  VA. 

i)  F.  M.  Guess  gave  an  invited  talk  on  "Mean  residual  life:  theory  and 
applications"  at  the  Southeastern  Regional  American  Mathematical 
Society  meeting,  Charlotte,  NC,  Octo^r  18,  1986. 

v)  K.  F.  Yu  and  W.  J.  Padgett  attended  Joint  Statistical  Meetings, 

Chicago,  August  17-21,  1986,  and  Dr.  Yu  presented  a  paper  entitled  "A 
note  on  a  renewal  theorem  for  a  moving  average  process." 

v)  F.  M.  Guess  emd  K.  F.  Yu  attended  the  9th  Annual  Meeting  of  the  S.C. 
Chapter  of  ASA.  April  10,  1987,  in  Col’inbia,  SC,  and  were  judges  of 
the  student  presentations. 


INVENTIONS,  PATENT  DISCLOSURES,  AND  APPLICATIONS  STEMMING  FROM  TOE  RESEARCH 
PROJECT 

No  inventions  or  patents  have  stemmed  from  this  research. 

The  results  reported  in  Section  4  have  wide  application  in  the  estimation 
and  assessment  of  reliedsility  and  maintenance  of  military  equipment.  The 
various  nonparametric  procedures  with  censored  data  developed  in  this  project 
allow  accurate  estimation  of  medieui  and  mean  lifetime,  percentiles,  densities, 
survival  probability,  and  other  lifetime  characteristics  without  assuming 
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particular  forms  of  the  life  distribution.  In  particular,  the  kernel-type 
quantile  estimator  will  provide  more  accurate  estimates  of  the  medieui  lifetime 
and  smaller  percentiles  of  a  piece  of  equipment  them  the  usual  product-limit 
estimators. 

8.  OTHER  PROFESSIONAL  ACTIVITIES 

During  this  reporting  period,  June  1,  1884,  to  May  31,  1987,  the 
investigators  have  been  involved  in  numerous  other  professional  activities  that 
are  intimately  related  to  the  research  efforts  on  this  grant.  Professor 
Padgett  has  refereed  27  manuscripts  for  Technometrics,  the  Journal  of  the 
American  Statistical  Association,  Communications  in  Statistics,  Statistics  amd 
Probability  Letters,  Pakistan  Journal  of  Statistics,  Biometrika,  IEEE 
Transactions  on  Reliability,  Journal  of  Multivariate  Analysis,  Journal  of 
Statistical  Planning  and  Inference,  Operations  Research,  and  the  Journal  of  the 
University  of  Kuwait  (Science).  Professor  Yu  refereed  six  articles  for  the 
Annals  of  Statistics,  Sequential  Analysis,  Statistics  and  Probability  Letters, 
JASA,  amd  the  Journal  of  Statistical  Plamning  and  Inference.  Professor  Guess 
has  refereed  one  article  for  each  of  the  Journal  of  the  Americam  Statistical 
Association  and  the  Naval  Research  Logistics  Quarterly.  In  addition,  five 
papers  were  reviewed  for  the  Zentrablatt  fiir  Mathematik  and  21  were  reviewed 
for  the  Mathematical  Reviews.  Also,  four  research  proposals  were  reviewed  for 
the  APOSR,  four  for  the  ARO,  and  one  for  the  NSF. 

Professor  Padgett  was  a  member  of  the  International  Editorial  Board  of  the 
Ccxanunications  in  Statistics-Theory  and  Methods  euid  was  an  Associate  Editor  of 
Naval  Research  Logistics  Quarterly  and  the  Journal  of  Statistical  Conpjtation 
sni  Simulation.  In  addition,  he  was  the  Program  Co-Chair  for  the  SRCOS-ASA 
Summer  Research  Conference  in  Statistics  held  in  Mobile,  AL,  on  June  15-20, 


1986.  He  also  chaired  an  invited  paper  at  the  Spring  Statistics  Meetings  in 
Atlanta,  GA,  in  March,  1986.  Professor  Spurrier  was  a  member  of  the 
International  Editorial  Board  of  the  Communications  in  Statistics  -  Statistical 
Reviews  euid  was  an  Associate  Editor  of  Technometrics. 


